#include "Mesh.h"
#include "Element.h"
using namespace std;
int main()
{
    Element<2> *e =new  Quadrilateral_4_Element();
    vector<double> v;
    v.push_back(e->phi_xi(-1,-1,1));
    v.push_back(e->phi_xi(1,-1,2));
    v.push_back(e->phi_xi(1,1,3));
    v.push_back(e->phi_xi(-1,1,4));
    v.push_back(e->phi_xi(-0.5,-1,5));
    v.push_back(e->phi_xi(0,-1,6));
    v.push_back(e->phi_xi(0.5,-1,7));
    v.push_back(e->phi_xi(1,-0.5,8));
    v.push_back(e->phi_xi(1,0,9));
    v.push_back(e->phi_xi(1,0.5,10));
    v.push_back(e->phi_xi(0.5,1,11));
    v.push_back(e->phi_xi(0,1,12));
    v.push_back(e->phi_xi(-0.5,1,13));
    v.push_back(e->phi_xi(-1,0.5,14));
    v.push_back(e->phi_xi(-1,0,15));
    v.push_back(e->phi_xi(-1,-0.5, 16));
    v.push_back(e->phi_xi(-0.5,-0.5,17));
    v.push_back(e->phi_xi(0,-0.5,18));
    v.push_back(e->phi_xi(0.5,-0.5,19));
    v.push_back(e->phi_xi(0.5,0,20));
    v.push_back(e->phi_xi(0.5,0.5,21));
    v.push_back(e->phi_xi(0,0.5,22));
    v.push_back(e->phi_xi(-0.5,0.5,23));
    v.push_back(e->phi_xi(-0.5,0,24));
    v.push_back(e->phi_xi(0,0,25));
    for(int i =0 ; i < v.size(); i++)
    {
	cout << v[i] << endl;
    }
};
